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ABSTRACT 

Probing the magnetism of the upper solar chromosphere requires measuring and modeling the 
scattering polarization produced by anisotropic radiation pumping in UV spectral lines. Here we 
apply PORTA (a novel radiative transfer code) to investigate the hydrogen Lyct line in a 3D model of 
the solar atmosphere resulting from a state of the art MHD simulation. At full spatial resolution the 
linear polarization signals are very significant all over the solar disk, with a large fraction of the field 
of view showing line-center amplitudes well above the 1 % level. Via the Hanle effect the line-center 
polarization signals are sensitive to the magnetic field of the model’s transition region, even when its 
mean field strength is only 15 G. The breaking of the axial symmetry of the radiation field produces 
significant forward-scattering polarization in Lya, without the need of an inclined magnetic field. 
Interestingly, the Hanle effect tends to decrease such forward-scattering polarization signals in most 
of the points of the field of view. When the spatial resolution is degraded, the line-center polarization of 
Ly a drops below the 1 % level, reaching values similar to those previously found in ID semi-empirical 
models (i.e., up to about 0.5%). The center to limb variation of the spatially-averaged polarization 
signals is qualitatively similar to that found in ID models, with the largest line-center amplitudes 
at /i = cosO ^ 0.4 {0 being the heliocentric angle). These results are important, both for designing 
the needed space-based instrumentation and for a reliable interpretation of future observations of the 
Lya polarization. 

Subject headings: polarization — scattering — radiative transfer — Sun: chromosphere — Sun: tran¬ 
sition region — Sun: surface magnetism 


1. INTRODUCTION 

Spectrographs on board sounding rockets (e.g., 
Roussel-Dupre 1982) and space telescopes (e.g., Curdt et 
al. 2008) have successfully measured the intensity profile 
(Stokes /(A)) of the hydrogen Ly a line radiation at var¬ 
ious positions on the solar disk, showing its strong emis¬ 
sion character. The center of the hydrogen Ly a emission 
profile originates at the base of the chromosphere-corona 
transition region (TR), where the kinetic temperature 
suddenly jumps from less than 10^ K to more than 10^ K 
and the plasma changes from partially to practically fully 
ionized. However, nobody has ever measured the linear 
polarization profiles (Q(A) and U{X)) of the on-disk Ly a 
radiation. Recent radiative transfer investigations using 
one-dimensional (ID) and two-dimensional (2D) mod¬ 
els of the solar atmosphere suggest that the absorption 
and scattering of anisotropic radiation in the upper solar 
chromosphere produce measurable linear polarization in 
the hydrogen Ly a line radiation of the solar disk, in the 
line core and in the line wings (Trujillo Bueno et al. 2011; 
Belluzzi et al. 2012; Stepan et al. 2012). Moreover, via 
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the Hanle effect the line-core polarization should react to 
the presence of magnetic fields in the solar TR, with sen¬ 
sitivity to field strengths between 10 and 100 G (Trujillo 
Bueno et al. 2011, 2012). 

The above-mentioned theoretical investigations on 
scattering polarization and the Hanle effect in the Ly a 
line of the solar-disk radiation have motivated the devel¬ 
opment of the Chromospheric Ly a Spectropolarimeter 
(CLASP), a sounding rocket experiment presently under 
development (Kobayashi et al. 2012; Kano et al. 2012; 
Kubo et al. 2014). The top priority of CLASP will be to 
measure, for the first time, the linear polarization pat¬ 
tern (i.e., (3(A)//(A) and U(A)//(A); hereafter, Q/I and 
U /1) within a lA wavelength interval around the Ly a 
line center. The second objective will be to use such ob¬ 
servables to constrain the magnetic field of the upper so¬ 
lar chromosphere, by modeling the modification that the 
Hanle effect is expected to produce on the line-center lin¬ 
ear polarization. The wing signals of the Q/I profile can 
only be modeled taking into account partial frequency 
redistribution (PRD) and quantum interference between 
the sublevels pertaining to the two upper J-levels of the 
Lya line (Belluzzi et al. 2012). Fortunately, as pointed 
out by Belluzzi et al. (2012), the line-core signals of the 
Q/I and U/I profiles can be approximately modelled ne¬ 
glecting J-state interference and frequency correlations 
between the incoming and outgoing photons in the scat¬ 
tering events (i.e., within the limit of complete frequency 
redistribution, or CRD). This is important because the 
Hanle effect in Ly a operates only in the line core. 

The solar atmosphere is highly inhomogeneous and dy¬ 
namic. Plane-parallel, horizontally homogeneous at mo- 
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Fig. 1.— Top left panel: Horizontal variation of the height in the 3D model atmosphere where the Ly a line-center optical depth along the 
vertical direction is unity (t^=i = 1). Top right panel: Strength of the model’s magnetic field at the atmospheric heights where t^=i = 1; 
the resulting average magnetic field strength is (B) ^ 15 G. Bottom panels: Inclination (left) and azimuth (right) of the magnetic field at 
the heights where t^=i = 1. Note that between the two clusters of opposite polarity the magnetic field has a loop-like structure with a 
dominant azimuth and a strength varying between a few gauss and 60 gauss. 


spheric models can be used to estimate the significance 
of the scattering line polarization signals and their sen¬ 
sitivity to the Hanle effect (Trujillo Bueno et al. 2011, 
2012). However, for making progress in our understand¬ 
ing of the magnetic and thermal structure of the outer so¬ 
lar atmosphere it is important to confront spectropolari- 
metric observations with spectral synthesis calculations 
in increasingly realistic three-dimensional (3D) models 
resulting from state-of-the-art radiation magnetohydro¬ 
dynamic (MHD) simulations, such as those carried out 
with the Bifrost code (Gudiksen et al. 2011). To make 
feasible such forward modeling research, Stepan & Tru¬ 
jillo Bueno (2013) developed PORTA, a 3D multilevel 
radiative transfer code for modeling the intensity and po¬ 
larization of spectral lines using massively parallel com¬ 
puters. 

Stepan et al. (2012) investigated the linear polariza¬ 
tion of the hydrogen Ly a line produced by scattering 
processes in a two-dimensional (2D) model of the solar 


atmosphere. This 2D model was a slice from the Bifrost 
3D model atmosphere used by Leenaarts et al. (2012) to 
investigate the intensity profiles of the H a line. Such 3D 
model of the extended solar atmosphere is representative 
of an enhanced network region and the aim of the present 
paper is to show in detail the results we have obtained by 
solving, with PORTA, the full 3D problem of the gener¬ 
ation and transfer of scattering polarization in the Ly a 
line, taking into account the Hanle effect caused by the 
model’s magnetic field. We do not show the polarization 
of the Zeeman effect because it produces rather negligi¬ 
ble signals in Ly a (e.g., Trujillo Bueno 2014). To study 
the full 3D problem of scattering polarization and the 
Hanle effect in chromospheric and transition region lines 
is important because 3D models of the solar atmosphere 
lack the translational symmetry of the previously inves¬ 
tigated ID and 2D cases, which provided only a glimpse 
on the complexity of the expected scattering polarization 
signals in the real 3D world. 
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Fig. 2. — Grotrian diagram of the hydrogen model atom. The 
components of Lya are indicated along with their Einstein 
coefficients (in s“^). The grey arrows connecting the fine-structure 
levels of n = 2 indicate the weakly inelastic depolarizing transitions 
due to proton and electron collisions. 

This paper is organized as follows. After formulat¬ 
ing the radiative transfer problem in Section 2, we show 
the radiation field tensors (Section 3) and the multipo¬ 
lar components of the atomic density matrix (Section 4) 
that result from the self-consistent solution in the cho¬ 
sen 3D model atmosphere. Section 5 considers the spa¬ 
tially resolved case, showing the fractional polarization 
of the emergent Ly a radiation for two scattering geome¬ 
tries (disk center and close to the limb) in the absence 
and in the presence of the model’s magnetic field. In 
Section 6 we study the spatially averaged polarization 
signals, paying particular attention to their sensitivity 
to the spatial resolution element, to their center to limb 
variation, to their magnetic sensitivity via the Hanle ef¬ 
fect, and to their detectability with a spectropolarimeter 
like CLASP. Section 7 summarizes our main conclusions. 

2. THE RADIATIVE TRANSEER PROBLEM 

The 3D model atmosphere considered in this paper is a 
snapshot of a radiation MHD simulation performed with 
the Bifrost code (Gudiksen et al. 2011) taking into ac¬ 
count non-equilibrium hydrogen ionization (see details 
in Carlsson et al. 2014). We used snapshot 385 of the 
simulation “en024048-hion”, which is the same simula¬ 
tion recently made publicly available through the In¬ 
terface Region Imaging Spectrograph project (IRIS; De 
Pontieu et al. 2014) at the Hinode Science Data Center 
Europe (http://sdc.uio.no; see IRIS Technical Note 33). 
This 3D model atmosphere, which has horizontal peri¬ 
odic boundary conditions, encompasses the upper part 
of the convection zone, the photosphere, chromosphere, 
transition region and corona, and is identical to that used 
by Leenaarts et al. (2012) to investigate the Ha intensity 
profile. The model’s magnetic field configuration is repre¬ 
sentative of an enhanced network region; it has magnetic 
field lines connecting two clusters of photospheric mag¬ 
netic concentrations with two dominant opposite polar¬ 
ity regions 8 Mm apart, which reach chromospheric and 
coronal heights (see Eig. 1). 

The Ly a line results from two blended transitions be¬ 
tween the ground level of hydrogen, ls‘^Si/ 2 , arid fhe 
and 2p‘^P^j2 excited levels. The Ly a scatter¬ 
ing polarization calculations of this paper were carried 
out fixing the model’s neutral hydrogen number density 
at each spatial grid point (which were computed with 
the Bifrost code taking into account time-dependent hy¬ 
drogen ionization) and solving the problem of the gen¬ 
eration and transfer of polarized radiation assuming sta¬ 
tistical equilibrium in a hydrogen model atom with only 


bound-bound transitions. The Bifrost ionization balance 
for hydrogen was computed using the method described 
in Leenaarts et al. (2007), which assumes that the Lyman 
transitions are in detailed balance. 

The hydrogen model atom we have used in the scat¬ 
tering polarization calculations (see Eig. 2) is composed 
of the above-mentioned levels and the 2s ‘^Si /2 level. The 
collisional coupling of this level with the P-levels may 
contribute to a slight depolarization of the Ly a line in 
regions where the density of perturbers (protons and elec¬ 
trons) is significantly larger than 10^^ cm“^ (Stepan & 
Trujillo Bueno 2011). Since in models of the upper solar 
chromosphere the proton and electron density is typi¬ 
cally smaller than lO^^crn”^, even a three-level model 
atom (i.e., without the 2s‘^Sii2 level of Eig. 2) could be 
sufficiently realistic for computing the line-center polar¬ 
ization of Lya. The only level that produces scattering 
polarization in the core of the Ly a line is the 2p ^Ps /2 up¬ 
per level. It is interesting to note that the suitability of 
the above-mentioned three-level model atom implies that 
the atomic polarization of the 2 p^P^j 2 level (population 
imbalances and quantum coherence between its magnetic 
sublevels) is dominated by the absorption of anisotropic 
radiation in the resonance line transition itself and by 
the action of the Hanle effect on such upper level, with¬ 
out any serious impact of collisional depolarization. The 
critical Hanle field of the level 2p ^i 3/2 is Bh^^^ G and 
the scattering polarization of the on-disk Ly a radiation 
is sensitive to magnetic strengths between 10 and 100 G, 
approximately (Trujillo Bueno et al. 2011) 

To compute the emergent Stokes profiles /, Q and U 
at each surface point of the atmospheric model we have 
to solve the following radiative transfer equations for any 
desired line of sight {X being /, Q or U) 

fx = Sx-X, ( 1 ) 

where dr = rjjds defines the monochromatic optical path 
length along the ray direction (with s the corresponding 
geometrical distance), r]i is the absorption coefficient, 
and Sx = ^xlvi is the source function of the Stokes 
parameter X, with ex the emission coefficient.^ While 
the emission coefficient e/ has contributions from the two 
blended transitions indicated in Eig. 2, eg and eu depend 
only on the atomic polarization of the 2p ^P ^/2 upper level 
(because it is the only Ly a level that can carry atomic 
alignment). We quantify the population and atomic po¬ 
larization of this j = 3/2 level by means of the multipo¬ 
lar components of its atomic density matrix (i.e., the pQ 
components, with X = 0, 2 and Q = —X,..., X, as quan¬ 
tified in the reference system whose quantization axis is 
along the local vertical) or, equivalently, by the following 
quantities (with the pQ components normalized to the 
fraction of neutral hydrogen atoms that are in the ground 
level): Sq = • lu our radiative transfer 

calculations with PORTA we have used the exact S'/, Sq 

^ Note that in our radiative transfer problem 77 / is the only 
non-zero coefficient of the propagation matrix of the Stokes-vector 
transfer equation, because the Zeeman splitting is negligible com¬ 
pared to the width of the hydrogen Ly ol line and the angular mo¬ 
mentum of its lower level is j = 1/2 and, therefore, cannot carry 
atomic alignment (i.e., we do not have “zero-field” dichroism here; 
cf. Trujillo Bueno & Landi Degl’Innocenti 1997). 
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Fig. 3.— Physical quantities across the vertical 2D slice a,t X = 12 Mm in the 3D model. Upper left panel: Kinetic temperature. The 
white solid line indicates the height, associated to each horizontal Y-point, where the Lya line-center optical depth is unity along a LOS 
with II = 1. Other panels: Components of the radiation field tensors normalized to Jq . 


and Su expressions, which take into account that the two 
Ly a transitions are not perfectly overlapped. Such tran¬ 
sitions are actually very close in wavelength (they are 
only about 6 niA apart) and, for simplicity, we give here 
the Sq and Su expressions for the perfect overlapping 
case: 

— A/3/i\/l — /i^(cos xRe[5'^] — sinxIm[S'^]) 

- ^(1-^M^)(cos2xRe[S'|] - sin2xIm[S'|])| , 

and 

Su= r — /i2(sinxRe[5'i] + cosxIm[S'f]) (3) 

-h/u(sin2xRe[S'|] -h cos2xIm[S'|])| , 

where 0 = arccos /i and x are the inclination with respect 
to the solar local vertical and azimuth of the ray, respec¬ 
tively, and r = r]\(j)y/irii = ri\(j)j ^/+ r]j) (with 77 J 
and ? 7 j the continuum and the line-integrated opacities, 
respectively, and the Voigt line absorption profile). 
In such equations the reference direction for Stokes Q is 
in the plane formed by the ray’s propagation direction 
and the vertical Z-axis. We point out that Eqs. ( 2 ) and 
(3), with their corresponding Eqs. (1), imply that the 


sources of Stokes Q are the population imbalances {Sq) 
and the quantum coherence (5'g, with Qy^O), while the 
Stokes U signal is produced only by quantum coherence 
(5'g, with Qy^O) between pairs of sublevels pertaining to 
the 279 ^/ 3/2 excited level. 

To solve Eqs. ( 1 ) for any desired line of sight (LOS) 
we need to know first 77 / and the Sx source function 
components at each spatial grid point. To obtain the 
self-consistent values of such quantities we have to solve 
jointly the statistical equilibrium equations for the mul¬ 
tipolar components (pg (j)) of the atomic density matrix 
corresponding to each atomic level j and the radiative 
transfer equations for each of the radiative transitions 
in the chosen model atom. The general expressions of 
such equations can be found in Sections 7.2 and 7.13.C of 
Landi Degl’Innocenti & Landolfi (2004), which we solved 
applying PORTA. We used a Gaussian quadrature for the 
ray inclinations with respect to the vertical Z-axis (with 
5 inclinations per octant) and the trapezoidal rule for the 
ray azimuth (with 4 azimuths per octant). In our radia¬ 
tive transfer calculations we took into account the im¬ 
pact of the Doppler shifts, produced by the macroscopic 
velocity gradients of the 3D atmospheric model, on the 
anisotropic radiation pumping. Although this physical 
ingredient is not very important for modeling the scatter¬ 
ing polarization of the hydrogen Ly a line (because the 
line is very broad), PORTA takes it fully into account 
because the dynamic state of the solar atmosphere may 
produce observable signatures on the scattering polar- 
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Fig. 4. — The fractional atomic alignment components, pg/pQ, 
of the 2p^P^j2 upper level, in the same vertical slice as in Fig. 3. 


ization of some chromospheric lines (Carlin et al. 2013). 
Our calculations with PORTA were carried out in the 
MareNostrum III supercomputer of the Barcelona Su¬ 
percomputing Center. The total computing time of the 
non-LTE calculations carried out for the investigation of 
this paper exceeded one million CPU hours. 

To grasp better the physics of this complex radiative 
transfer problem it is useful to write down a simplified set 
of equations, instead of the general equations for the four- 


level model atom of Fig. 2 we have actually solved with 
PORTA to obtain all the results shown in this paper. 
To this end, we note that the weakly inelastic collisions 
with protons and electrons that couple the 2s ^5^/2 level 
with the 2 p‘^Pij 2 and 2 p‘^P^j 2 ones, which we have taken 
into account in our 3D radiative transfer calculations, do 
not really have a significant depolarizing effect on the 
Lya polarization in the considered model atmosphere 
(Stepan & Trujillo Bueno 2011). Taking also into ac¬ 
count that elastic collisions with neutral hydrogen atoms 
can be safely neglected when modeling the line-core po¬ 
larization of Ly a and that the only level contributing to 
the Lya scattering polarization is the 2 p‘^P^j 2 one, we 
have the following approximate expressions for the Sq 
tensors that quantify the excitation state of the 2 p^P^j 2 
upper level (Manso Sainz & Trujillo Bueno 2011): 


5o° = (l-e)Jo° + eB, 


( 4 ) 


(Si \ 




Re[S'i2] 


Re[J?] 


lm[Sf] 

(1-e) 

V2 



Re[5|] 

Re[Ji] 







( 

0 M\2 Afi3 0 


+ Tm (1 “ e) 


M21 0 M23 M24 Af 25 

M32 0 M34 M35 

0 M 42 Af43 0 Af45 

\ 0 M52 M53 M54 0 y 


( 5 ) 

/Si \ 

Re[5i2] 

Re[5|] 

\lm[Sl]) 


where is the Planck function at the By a line-center 
frequency, e is the photon collisional destruction proba- 
blility, Tu = 1.869 x (with the magnetic strength 

B in gauss), and the Mij coefficients of the magnetic ker¬ 
nel M depend on the inclination (Ob) of the local mag¬ 
netic field vector with respect to the vertical Z-axis and 
on its azimuth {xb)- The frequency averaged (over the 
absorption profile (j)^,) radiation field tensors Jq that ap¬ 
pear in such equations quantify the mean intensity (Jq)? 
the anisotropy {Jq)^ and the breaking of the axial sym¬ 
metry (the real and imaginary parts of and ) of the 
incident radiation field at each point within the medium; 
their explicit expressions can be found in Manso Sainz 
& Trujillo Bueno (2011), but note that the Stokes I of 
the incident Lya radiation has to be calculated taking 
into account that it results from two blended transitions 
(see Fig. 2). The equations for the Sq (atomic) quan¬ 
tities have a clear physical meaning. In particular, note 
that the magnetic operator M couples locally the K = 2 
components among them (the Hanle effect). 


3. THE RADIATION FIELD TENSORS 

The top left panel of Fig. 3 shows the spatial variation 
of the kinetic temperature across a 2D slice through the 
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Fig. 5. — Synthetic disk center observation (LOS with = 1). Top left panel: the calculated emergent intensity at the Lya line-center, 
quantified in terms of a radiation temperature. Top right panel: the total fractional linear polarization P = ■sjQ‘^ + U‘^ jI at the Lya line 
center. Bottom panels: the corresponding Q/I and U/I line-center signals, with the positive reference direction for Stokes Q along the 
T-axis. 


above-mentioned 3D snapshot model atmosphere; it cor¬ 
responds to the Z-Y plane at X = 12 Mm in Fig.l. Note 
that the height of the TR changes drastically from one 
horizontal location to another, outlining a highly crum¬ 
pled surface. As shown by the white solid curve, the 
corrugated surface that delineates the model’s TR prac¬ 
tically coincides with that defined by the Ly a line-center 
optical depth unity along the line of sight (LOS). 

The other panels of Fig. 3 show the radiation field 
tensors Jq of the Lya radiation, normalized to Jq. 
Such quantities quantify the fractional anisotropy and 
the symmetry properties of the incident radiation at each 
point within the medium. We calculated them after ob¬ 
taining with the 3D radiative transfer code PORTA the 
self-consistent solution in the 3D snapshot model illus¬ 
trated in Fig. 1, using the four-level model atom of Fig. 2 
and taking into account the (very weak) collisional cou¬ 
pling between the 2s‘^Sii2 fFe two P levels. Of 
particular interest is the spatial variation of the frac¬ 
tional anisotropy (top right panel); it is negligible in the 
photosphere and chromosphere of the model, but it sud¬ 
denly becomes significant right at the location of the TR 
with negative values of about 10%. The other radia¬ 
tion field tensors shown in the remaining panels quantify 
the breaking of the axial symmetry of the Ly a radiation 
at each point within the model atmosphere; note that 
at the atmospheric heights where the line center optical 


depth is unity such tensors may have sizable positive or 
negative values, depending on the horizontal point under 
consideration. 

4. THE ATOMIC LEVEL POLARIZATION 

The multipolar components of the atomic density ma¬ 
trix (i.e., the pQ quantities or, equivalently, the Sq ones) 
quantify the excitation state of the atomic level under 
consideration. As mentioned above, the scattering po¬ 
larization of Ly a is produced by the atomic polarization 
of its 2 p‘^P^j 2 upper level; therefore, in this section we 
provide information on the spatial variation of the level’s 
fractional atomic polarization throughout the same Z-Y 
plane considered in Eig. 3. Eqs. (5) indicate that for 
Tu = 0 (i.e., in the unmagnetized case) each multipolar 
component pg of the 2p ^Ps /2 level is directly proportional 
to its corresponding radiation field tensor component. 
We can say that there is a transfer of ‘order’ from the ra¬ 
diation field to the atomic system, and Eqs. (l)-(3) show 
that the atoms of a polarized (‘ordered’) atomic system 
emit linearly polarized radiation. Eqs. (5) indicate that 
if all Jq were zero then all S'g would vanish, even in the 
presence of a magnetic field. Through the Hanle effect 
the magnetic field modifies the polarization of the atomic 
level according to the F^^ term of Eqs. (5), and this has 
an impact on the emergent linear polarization Q/I and 
U/I signals. 
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Fig. 6. — For the forward scattering case (LOS with ii = 1), the top panels show the total fractional linear polarization P = a/ -h U‘^ jI 
at the Lya line center assuming the absence of magnetic fields (top left panel) and considering the Hanle effect produced by the model’s 
magnetic field (top right panel). The bottom left panel shows the P image that results when assuming that at each point within the 3D 
model the magnetic strength is in the Lyo; line Hanle saturation regime (i.e., B > 250 G). The bottom right panel shows the difference 
between the true total fractional linear polarization (obtained by computing /, Q and U taking into account the Hanle effect of the model’s 
magnetic field) and that corresponding to the zero-field reference case. 


Fig. 4 shows the multipolar components of the density 
matrix corresponding to the 2 p‘^P^j 2 level, normalized to 
the po value of the same level (which is proportional to 
its overall population). Note that the Po/Po P^iiel (which 
quantifies the population imbalances among the sublevels 
of the 2 p^i 3/2 level) shows mainly negative values at the 
atmospheric heights where the Lya line center optical 
depth is unity along the LOS, while each of the quantum 
coherence shown in the other panels shows both positive 
and negative values. 

5. THE SPATIALLY RESOLVED Q/I AND U/I 
SIGNALS 

This section considers the case of spatially resolved 
observations. It shows maps of the line-center scatter¬ 
ing polarization, which we have obtained by calculat¬ 
ing the emergent Stokes profiles of the Lya radiation 
after obtaining the self-consistent solution of the above- 
mentioned 3D radiative transfer problem. Of particu¬ 
lar interest is to investigate the spatial variability of the 
emergent Q/I and U/I line-center signals and their sen¬ 
sitivity to the model’s magnetic field. To this end, we 
solved the 3D radiative transfer problem of resonance line 
polarization for the following two cases: (1) taking into 
account the Hanle effect produced by the magnetic field 


of the 3D snapshot model atmosphere and (2) assuming 
H = 0 G at each spatial grid point of the 3D model. We 
show the results for two line of sights (LOS) defined by 
p = cos6> (with 0 the heliocentric angle): p = 1 (corre¬ 
sponding to a disk center observation) and p = 0.4 (cor¬ 
responding to another on-disk observation, but at about 
80 seconds of arc from the solar limb). In all our cal¬ 
culations of the emergent Stokes profiles the azimuth of 
the corresponding LOS was X = 0° - that is, the LOS 
was contained in the X-Z plane. In all the following fig¬ 
ures we have chosen the Y-axis as the positive reference 
direction for Stokes Q, and for any LOS with p < 1 we 
chose this Y-direction as the parallel to the nearest limb. 

Eig. 5 shows, for the forward scattering case of a 
disk center observation (LOS with p = 1), the inten¬ 
sity (top left panel), the total fractional linear polariza¬ 
tion (P = ^/Q^X-TP/I, top right panel), and the Q/I 
(bottom left panel) and U/I (bottom right panel) sig¬ 
nals at the center of the emergent Lya intensity pro¬ 
file. The calculated Lya line-center intensity variations 
across the surface of the 3D atmospheric model show 
fibril-like structures, reminiscent of those observed at 
1 arcsecond resolution by the VAULT sounding rocket 
(e.g., Vourlidas et al. 2010). Although we do not show 
illustrations here, we point out that the effects of 3D 
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Inclination of the normal vector of the transition region 
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Fig. 7. — The inclination On of the vector that is normal to 
the corrugated surface that delineates the model’s transition re¬ 
gion, defined here as that where the Ly a line-center optical depth 
along the vertical direction is unity (Ty^=i = 1). The dark areas 
correspond to the places where the model’s transition region is 
horizontal, as in a ID model atmosphere, whereas the bright areas 
indicate the places where the TR is vertical. The average inclina¬ 
tion is (^n) = 26°. We note that a very similar figure is obtained 
if the normal vectors to the TR are defined as being parallel to the 
gradient of the kinetic temperature at the line-center formation 
heights. 


radiative transfer do have an impact on the emergent 
line-center intensities, as expected from the smoothing ef¬ 
fect produced by horizontal radiative transfer under non- 
LTE conditions (e.g., Kneer 1981); the so-called 1.5D ap¬ 
proximation gives a larger line-center intensity contrast 
across the field of view^. The same applies to other spec¬ 
tral lines such as hydrogen Ha for which this effect is 
even more significant due to the lower opacity of the line 
(Leenaarts et al. 2007). 

It is interesting to note that at high spatial resolution 
the linear polarization signals of the emergent Ly a line- 
center radiation are very significant (see Fig. 5). At first 
sight this result might appear surprising, given that in 
Fig. 5 we are considering the forward scattering case of a 
disk-center observation. However, in this scattering ge¬ 
ometry zero linear polarization should be expected only 
for the idealized case of an unmagnetized and horizon¬ 
tally homogeneous stellar atmosphere model. As soon 
as there is a magnetic field inclined with respect to the 
symmetry axis of the incident radiation field (hereafter, 
inclined field) and/or horizontal atmospheric inhomo¬ 
geneities (e.g., in the gas kinetic temperature and/or den¬ 
sity) the ensuing symmetry breaking effects can produce 
significant line scattering polarization signals (see Manso 
Sainz & Trujillo Bueno 2011; Stepan et al. 2012, and 
more references therein). The linear polarization signals 
of Fig. 5 are due to the breaking of the axial symmetry 
of the Ly a radiation at each spatial point within the 3D 
model atmosphere, as quantified by the real and imag¬ 
inary parts of the and radiation field tensors, as 
well to the magnetic field of the atmospheric model (see 

® With this approximation, every column of the model atmo¬ 
sphere is approximated by a ID plane-parallel model, so that hor¬ 
izontal radiative transfer effects are neglected. 
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Fig. 8. — The same as in Fig. 5, but for the case of a close to the 
limb observation (p, = 0.4). The visualization takes into account 
the foreshortening effect corresponding to the chosen inclined line 
of sight. 


Eqs. 4-5). The Q/I and U/I line-center signals shown in 
the bottom panels of Fig. 5 have similar values that fluc¬ 
tuate in sign across the field of view, as expected from 
the particularization of Eqs. (2) and (3) to the case of a 
LOS with /i = 1. In particular, the image corresponding 
to the total linear polarization at the line center (see the 
top right panel of Fig. 5) shows an impressive detail of 
fine-scale structuring. 

As explained above, the forward-scattering polariza¬ 
tion signals are caused by both the model’s horizon¬ 
tal atmospheric inhomogeneities in temperature, den¬ 
sity, and macroscopic velocity, and by the magnetic field. 
As shown in the top left panel of Fig. 6, there are 
very significant forward scattering polarization signals 
in the absence of magnetic fields. Such zero-field lin¬ 
ear polarization signals are solely due to the symmetry 
breaking caused by the horizontal atmospheric inhomo- 
geneties; therefore, their very existence indicates that the 
so-called 1.5D approximation is unsuitable for modeling 
high spatio-temporal resolution observations of the scat¬ 
tering polarization in Lya. Fig. 6 shows that there is 
a significant difference between the true total fractional 
linear polarization (calculated taking into account the 
Hanle effect produced by the magnetic field of the 3D 
model atmosphere) and that corresponding to the zero- 
field reference case. As expected, the action of the Hanle 
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Fig. 9. — Histogram of the total fractional linear polarization 
P = + U‘^/1 at the Lyo; line center. Top panel: Disk center 

case (LOS with = 1). Bottom panel: close to the limb case (LOS 
with II = 0.4). 

effect in forward scattering geometry is found mainly 
along the diagonal of the image seen in Fig. 6, where 
the magnetization of the model atmosphere is the largest. 
However, as seen in the figure, the model’s magnetic field 
tends to depolarize the emergent radiation, and this oc¬ 
curs mainly in the strongly polarized regions of the un¬ 
magnetized case. This contrasts with the case of a plane- 
parallel model atmosphere, where in the absence of mag¬ 
netic fields the symmetry axis of the incident radiation 
coincides with the vertical direction and the Hanle effect 
of an inclined magnetic field creates linear polarization in 
forward scattering geometry but depolarization in close 
to the limb observations. In our 3D model atmosphere 
the symmetry axis of the radiation field that illuminates 
each of the points of the TR does not in general coincide 
with the local vertical. In Fig. 7, we attempt to visualize 
the geometrical complexity of the present 3D problem by 
showing, at each point of the field of view, the inclina¬ 
tion of the vector that is normal to the TR surface. In 
the unmagnetized case the strongly polarized regions are 
typically those for which the inclination of the local nor¬ 
mal vector to the TR surface is larger, and depolarization 
by the Hanle effect is statistically more likely. In regions 
where the local normal vector to the TR surface is pre¬ 
dominantly vertical the forward scattering signals of the 
unmagnetized case are typically lower and the Hanle ef¬ 
fect produced by the model’s magnetic field tends to be 
more similar to that found in the case of a plane-parallel 
model atmosphere. These two situations can be under¬ 


stood by analogy with the scattering polarization signals 
calculated in a ID plane-parallel model near the limb and 
at the disk center, respectively (cf., Trujillo Bueno et al. 
2011 ). 

We now consider a LOS with /r = 0.4, which corre¬ 
sponds to the case of an observation relatively close to the 
solar limb. Fig. 8 shows maps of the line-center intensity 
(top left panel), of the total fractional linear polariza¬ 
tion (top right panel), as well as of the Q/I (bottom left 
panel) and U/I (bottom right panel) line-center signals. 
The particularization of Eqs. (2) and (3) to the LOS un¬ 
der consideration (or to any LOS with /i<l) indicates 
that eu (the emissivity in Stokes U at each point along 
the LOS) depends only on the quantum coherence {p\ 
and P 2 ) of the 2 p‘^P^j 2 excited level, while eg (the emis¬ 
sivity in Stokes Q at each point along the LOS) depends 
also on the population imbalances between the sublevels 
of the excited level. As indicated by Eq. (5), the popula¬ 
tion imbalances, quantified by Sq, are determined mainly 
by the radiation field tensor Jq, which at the crumpled 
surface that delineates the model’s TR tends to be neg¬ 
ative (as in the ID atmospheric models considered by 
Trujillo Bueno et al. 2011). We therefore expect that the 
Q/I signals of the emergent Ly a radiation across the 
surface of the 3D model have mainly negative sign. Eor 
the U/I signals we expect the occurrence of positive and 
negative values to be equally likely, following the fact 
(see Eq. 5) that the radiation field tensors Jf and 
(which quantify the symmetry breaking of the incident 
radiation field) play a key role on the quantum coherence 
quantified by and 

The upper panel of Eig. 9 shows the histogram of the 
line center polarization amplitudes of the emergent Lya 
radiation calculated in forward scattering geometry (LOS 
with /i = 1) assuming B = OG. Most of the linear po¬ 
larization signals have values between 0 and 1 %, but 
there is a clear tail of stronger polarization signals ex¬ 
tending up to about 3%. In this forward scattering ge¬ 
ometry the mean fractional linear polarization amplitude 
at the maximum spatial resolution of the 3D model is 
«0.7% 

The bottom panel of Eig. 9 shows an analogous his¬ 
togram, but for the case of a LOS with p = 0.4. As in 
the /i = 1 case (see the upper panel of Eig. 9), almost all 
the surface points are significantly polarized. In contrast, 
in this close to the limb scattering geometry there is only 
a minority of surface points with P < 1 % and there is a 
significantly larger fraction of points of the field of view 
showing strong polarization signals with P > 3%. As a 
result, the mean fractional linear polarization amplitude 
at the maximum spatial resolution of the 3D model is 
now P^=o .4 ~ 1.5%. 

Eig. 10 shows, for the case of a close to the limb LOS 
with ji = 0.4, the synthetic Q/I and U/I profiles com¬ 
puted at each T-point without accounting for any spa¬ 
tial degradation; the T-direction here can be assumed 
to be along the spectrograph’s slit, assuming it to be 
oriented parallel to the nearest limb (i.e., the Y axis). 
Of great interest is to compare the top panels (the zero- 
field reference case) with the bottom panels (where the 
Hanle effect of the model’s magnetic field has been taken 
into account). In both cases the Q/I signals are mainly 
negative, as expected from the behavior of the radiation 
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Fig. 10. — Simulated close to the limb observation {ji = 0.4). Emergent Q/I and U/I profiles calculated at each point along the E-axis 
direction at X=12 Mm in Fig. 1. They have been obtained by solving the full 3D radiative transfer problem, neglecting (top panels) and 
taking into account (bottom panels) the Hanle effect produced by the model’s magnetic field. The positive reference direction for Stokes 
Q is parallel to the nearest limb (i.e., along the E-axis). 



Fig. 11. — Average total fractional polarization signal P for 
the forward-scattering case of a disk-center observation as a func¬ 
tion of the spatial resolution R (in arcseconds) of a square R x R 
element. The dotted lines indicate the resolution for which the 
average polarization signal is reduced by a factor 1/2 with respect 
to the full-resolution case. 


field anisotropy {Jq/Jq) illustrated in Fig. 3. However, 


the polarization amplitudes are larger in the absence of 
magnetic fields, with a sizable factor of difference in the 
atmospheric regions where the magnetic field of the 3D 
model is relatively strong and inclined (e.g., see the re¬ 
gion around Y = 15 Mm). The action of the Hanle effect 
can also be clearly seen in the U/I signals (right panels). 

The results shown above indicate that at high spatial 
resolution we should expect conspicuous linear polariza¬ 
tion signals in the Ly a line core, with amplitudes often 
larger than 1%. These Q/I and U/I signals are pro¬ 
duced by anisotropic radiation pumping in the upper so¬ 
lar chromosphere and they are modified, via the Hanle 
effect, by the vector magnetic field of the model’s TR. In 
the following section, we turn our attention to the case 
of spatially averaged observations. 

6. THE SPATIALLY AVERAGED SIGNALS 

The solar atmosphere is a highly inhomogeneous and 
dynamic medium, much more complex than the idealiza¬ 
tion of a ID, static, plane-parallel semi-empirical model 
atmosphere. In a ID model atmosphere with axial sym¬ 
metry around the local vertical the only possibility to 
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Fig. 12. — The calculated CLV of the Qjl signal that results 
from spatially averaging the Stokes / and Q line-center signals at 
each 11 location, averaged over all the LOS azimuths x? with (solid 
curve) and without (dotted curve) the Hanle effect produced by 
the model’s magnetic field. The dashed curve corresponds to the 
Hanle effect saturation regime. The dash-dotted curve shows the 
calculated CLV in the semi-empirical model C of Fontenla et al. 
(1993). 



Fig. 13.— For a LOS with p, = 0.4 the figure shows the spatially 
averaged Q/I amplitude vs. the mean magnetic field strength of 
the model’s transition region. The arrow indicates the Q/I signal 
corresponding to the original 3D model. 

break the axial symmetry of the incident radiation field 
at each point within the medium is by the presence of 
a magnetic field inclined with respect to the solar ra¬ 
dius vector through the point under consideration (see 
Ishikawa et al. 2014, and more references therein). In the 
real solar atmosphere we have, in addition, the symme¬ 
try breaking produced by the horizontal inhomogeneities 
of the atmospheric plasma, in its density, velocity and 
kinetic temperature. These (non-magnetic) causes of 
symmetry breaking produce Q/I and U/I signals which 
could be confused with those resulting from the Hanle 
effect. 

In this paper we represent the complexity of the TR 
of the Sun by a single snapshot taken from the three- 
dimensional MHD simulation described in Section 2. One 
possible way to obtain quantitative information on the 
magnetization of the solar TR is to study the statisti¬ 
cal properties of the hya polarization signals, and their 
relation with the physical quantities that describe the 
3D atmospheric model. The first step towards this goal 
is to consider the average polarization signals and how 


they can be exploited to infer information on the aver¬ 
age magnetic field strength in the TR. In this section, we 
study the sensitivity of the emergent Q/I and U/I line- 
center signals to the model’s magnetic field, considering 
spatially averaged synthetic observations. As shown be¬ 
low, such spatially-averaged signals can be used to obtain 
quantitative information on the average magnetization of 
the solar TR. 

If not stated otherwise, the results of this section are 
shown for a LOS with azimuth X = 0°, i.e., with the Y- 
axis parallel to the nearest solar limb which is also our 
reference direction for Stokes Q. In order to simplify the 
notation, in what follows the spatially averaged quanti¬ 
ties {Q)/{I), {U)/{I), and (P) will be denoted by Q/I, 
U/I, and P. 

6.1. Spatial variability of the Lya signals 

If the spatial resolution of the observation is reduced, 
the emergent Stokes profiles /, Q, and U are averaged 
over the corresponding resolution element. At the disk 
center, the Q/I, U/I, and P signals are therefore gener¬ 
ally degraded because of cancellation of the positive and 
negative signals coming from different regions of the at¬ 
mosphere. Let us consider a square resolution element of 
size Rx R located at the solar disk center. In Fig. II we 
show the variation with R of the total polarization signal 
P at the Lya line center. This result has been obtained 
statistically, after computing the emergent Stokes profiles 
for many possible realizations of the resolution element 
position. As seen in the figure, at P = 3" the total polar¬ 
ization line-center signal decreases by a factor two with 
respect to the maximum resolution value. Three seconds 
of arc can therefore be regarded as a characteristic spatial 
scale of the polarized structures seen at the center of the 
Ly a line. We can therefore argue that unless the spatial 
resolution of the observation is significantly better than 
3 arcseconds, the polarization of the emergent radiation 
will be a mixture of signals coming from locations hav¬ 
ing significantly different plasma structures. Likewise, 
we can expect that for inferring the average properties 
of the solar TR, the spatial resolution of a single obser¬ 
vation should be much worse than 3 arcseconds, so that 
the influence of the local inhomogeneities is sufficiently 
reduced. 

For the case of a disk-center observation, the line- 
center signal P tends to zero in the limit of no spatial 
resolution. For any LOS with ja < 1, the spatially aver¬ 
aged signals will be generally non-zero and converging, in 
the limit of no spatial resolution, to the values given by 
the center-to-limb variation (CLV). This applies to the 
unmagnetized case, but as shown in Fig. 12 it happens 
also when we take into account the Hanle effect produced 
by the weak magnetic field of the chosen 3D model at¬ 
mosphere. 

6.2. Center-to-limb variation 

The CLV of the emergent spectral line radiation con¬ 
tains useful information on the structuring of the solar 
atmosphere. In quiet regions of the solar surface the 
observed CLV of the Stokes-/ profile of the hydrogen 
Lya line is known to be negligible (e.g., Roussel-Dupre 
1982). However, radiative transfer calculations in the 
semi-empirical model C of Fontenla et al. (1993) (here¬ 
after, FAL-C model), show that the CLV of the Q/I line- 
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center signal is rather significant (cf. the Fig. 2 of Tru¬ 
jillo Bueno et al. 2011). Here we study the CLV of the 
spatially-averaged Q/I line center signals in the chosen 
3D model. 

The results are shown in Fig. 12 for the case of az- 
imuthally averaged line of sights of inclination ji. Only 
the Q/I line-center signal is shown because, as expected, 
the spatially averaged U/I signal is practically zero for 
all heliocentric angles. Interestingly, the behavior of the 
resulting CLV of the line-center Q/I signal is qualita¬ 
tively similar to that found in the FAL-C model (see 
the dash-dotted line in Fig. 12). The largest fractional 
polarization amplitude is found around /i ~ 0.4 and it 
decreases towards larger heliocentric angles. In contrast 
to the CLV behavior found in ID models, the decrease 
of the Q/I line-center amplitude when approaching the 
limb of the 3D atmospheric model is not so significant 
and in the chosen 3D model the Q/I signal is still siz¬ 
able at near-limb locations. 

Without spatio-temporal resolution, the resulting Q/I 
signals drop below the 1% level, with amplitudes similar 
to those found in ID semi-empirical models. Therefore, 
their measurement requires advanced instrumentation, 
capable of reaching a polarimetric sensitivity of at least 
0.1%. Via the Hanle effect, the spatially-averaged Q/I 
signals are clearly sensitive to the model’s magnetic field 
(see Fig. 12). 

6.3. Magnetic sensitivity of the average spectrum 

Following the approach introduced in Section 4 of 
Stepan et al. (2012), we study the magnetic sensitivity 
of the spatially-averaged polarization amplitudes in the 
3D snapshot model atmosphere. We define the quantity 
(B) as the average magnetic field strength at the points 
where the optical depth at the Ly a line center is unity for 
the case of a disk-center observation. Given the “height 
of formation” of the Lya line center, (B) can also be 
understood as the average magnetic field strength in the 
TR. We follow the approach of Stepan et al. (2012) and 
multiply the magnetic field strength at each grid point 
of the 3D model by a dimensionless factor / > 0. In this 
way, we can study the dependence of the line polariza¬ 
tion on the average field strength f{B). The result for 
the Q/I line center amplitude is shown in Fig. 13 for a 
near-limb observation with fi = 0.4.^ As expected, the 
maximum polarization signal is found in the unmagne¬ 
tized case (/ = 0). The arrow in the figure indicates the 
actual value of the average TR magnetization (the / = 1 
case) which is (H) = 15 G. The Q/I polarization signal 
gradually decreases as the magnetic field increases. In 
contrast to the 2D result shown in Fig. 4 of Stepan et al. 
(2012), in the 3D model the Q/I signal is lower and its 
maximum amplitude is more similar to that found in the 
FAL-C unmagnetized model. 

6.4. Detectability of the signals 

The Chromospheric Ly a Spectropolarimeter (CLASP) 
has been designed for measuring the core of the Stokes 
/, Q and U profiles using a spectrograph with a slit of 
400", with a spatial resolution of about 3", a spectral res¬ 
olution of 0.1 A and a polarimetric sensitivity of 0.1% at 

^ The variation of U/I is not shown because this quantity re¬ 
mains practically zero for all values of /. 



0.1 1.0 10.0 100.0 
Averaging length [arcsec] 


Fig. 14. — Dependence of the slit-integrated Q/I amplitude 
on the averaging length along the spectrograph’s slit, with the slit 
located at /r = 0.4 and oriented parallel to the model’s solar limb. 
The figure shows the polarization signal modified by the Hanle 
effect of the model’s magnetic field (red solid curve), the solution 
without the magnetic field (dash-dotted curve), and the solution 
obtained with the magnetic field strength increased by a factor 4 
at every point of the atmospheric model (green solid curve). The 
error bars indicate the standard deviation of the Q/I amplitudes 
calculated from the statistics of randomly positioned slits. 



Fig. 15. — Amplitude of the (Q)/(I) signal, integrated along 
the 400" slit, as a function of the averaged magnetic field strength 
(B) = /{Bsd) of the model’s transition region. Blue curve: Slit 
oriented radially and extending from /r = 0.1 to = 0.8. Red 
curve: Slit oriented tangentially with the slit center located at 
II = 0.4. 

the 3(7 level (see Kobayashi et al. 2012; Kano et al. 2012; 
Kubo et al. 2014). In order to reach this level of po¬ 
larimetric sensitivity it will be necessary to significantly 
sacrifice the temporal resolution and to partially aver¬ 
age the observed Stokes /, Q, and U signals along the 
direction of the spectrograph’s slit, in order to increase 
the signal-to-noise ratio. Nevertheless, in addition to the 
important goal of measuring for the first time scattering 
polarization in a FUV line (hydrogen Lya) whose core 
originates in the solar TR, CLASP can still aim at con¬ 
straining the average magnetization of the TR by partial 
spectrum averaging along the spectrograph slit. 

An important question to consider is whether the spec¬ 
trograph’s slit, which selects solar light along a distance 
of 400” on the solar disk (i.e., almost one half of the 
apparent solar radius), should be oriented radially or 
tangentially with respect to the solar limb, in order to 
maximize the chances of detection of the expected scat- 
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tering polarization in Lya and the scientific output of 
the CLASP experiment. We dedicate this section to an¬ 
alyze the pros and cons of both options. Although both 
options seem to be suitable for detecting the scattering 
polarization signals, we have some reasons to prefer the 
radial option (see below). 

We begin our analysis by considering the case of tan¬ 
gential orientation of the spectrograph’s slit, with its cen¬ 
ter located at the heliocentric angle /i = 0.4, i.e., at an 
on-disk position where the spatially-averaged Q/I sig¬ 
nal is maximal (cf. Fig. 12). We have studied the de¬ 
pendence of the spatially-averaged Q/I line-center signal 
on the length L (in arc seconds) along the slit’s direc¬ 
tion over which we average the computed Stokes /, Q 
and U signals. We have done this analysis by averaging 
the individual signals corresponding to a large number 
of randomly positioned virtual slits on top of the 3D 
model grid. This has been done by selecting randomly 
the [X, F] point that corresponds to the first point of 
the slit. The average is then obtained by integration 
of the /, Q, and U signals along the randomly selected 
slit. If the virtual slit crosses the model’s boundary, then 
the integration continues along the same direction but 
starting at the opposite boundary, with a random po¬ 
sition along this boundary. This way, it is possible to 
obtain slit-integrated signals for a large number of slits 
that are longer than the model’s dimensions (which is 
about 33" X 33"). We have applied this procedure for 
averaging lengths between L = 0.2" and L = 400". The 
results are shown in Fig. 14 for three different degrees 
of the model’s magnetization. It is important to em¬ 
phasize that the spatialy-averaged Q/I signal depends 
on L: stronger signals are found, on average, for smaller 
lengths L. This is due to the fact that Q/I depends on 
correlations between Q and / at particular pixels. In our 
3D model, the regions with lower line intensity are, on 
average, more polarized than the high-intensity regions. 
Consequently, stronger signals are obtained for smaller 
lengths L. The spread of the resulting Q/I signals is in¬ 
dicated by the error bars in the figure, which correspond 
to the standard deviation of each Q/I value for any given 
length L. As expected, the larger L, the smaller the 
spread of the values. Averaging over the whole 400" dis¬ 
tance of the CLASP spectrograph’s slit would provide 
a 2cr test for detection of the zero magnetization versus 
the actual model magnetization {{B) = 15 G) of the 3D 
model. Similarly, L = 30" would be sufficient for a test 
distinguishing between a non-magnetized and a strongly 
magnetized {f{B) = GOG) transition region, with a sim¬ 
ilar level of confidence. 

Fig. 15 shows, both for the radial and tangential slit 
options, how the spatially-averaged (including averaging 
along the whole length of the spectrograph’s slit) Q/I 
line-center signal of the emergent Ly a radiation from the 
3D model varies with the average magnetization value of 
the model’s TR. For each average magnetization, f{B), 
a very large number of Q/I signals were obtained by 
a method analogous to the one described in the previ¬ 
ous paragraph. In the radial-slit option case, the fact 
that the Q/I signal varies with the heliocentric angle (see 
Fig. 12) has been taken into account. The amplitude of 
the Q/I signal predicted for the tangential configuration 
is slightly larger than in the radial case, but the difference 
is not very important and both slit orientations appear to 


be suitable for detecting the Ly a scattering polarization. 

Our reasons in favor of the radial slit option for a 
sounding rocket experiment such as CLASP are the fol¬ 
lowing. First, in the case of the tangential slit orienta¬ 
tion we would have data only for a small interval of fi 
values, for example with /i^O.4 at the central slit point 
(where the polarization amplitude is the largest). Obvi¬ 
ously, any significant pointing error would compromise 
the measurement. In contrast, the radial slit option pro¬ 
vides much more information on the CLV of the scat¬ 
tering polarization signals, and having such CLV data 
would facilitate reaching more reliable conclusions via 
detailed confrontations with radiative transfer modeling. 
Note that with a 400" radially-oriented slit one can have 
information from the solar limb till about fi = 0.8. With 
the radial slit option such CLV information would be 
available both at the line center (where the Hanle effect 
in Lya operates) and in the nearby wings of the Q/I 
profile (where effects of partial frequency redistribution 
and J-state interference come into play; see Belluzzi et al. 
(2012)). As shown by these authors, having in addition 
information on the CLV of the Q/I wing signals may help 
to constrain (together with the observed Stokes I pro¬ 
files) the thermal structure of the upper chromosphere. 
Another reason to prefer the radial slit option is related 
with the fact that with the CLASP (which will mea¬ 
sure the Lya radiation ±0.5 A from the line center) it 
will be difficult to determine with sufficient precision the 
zero offset of the linear polarization scale, and having the 
above-mentioned CLV information may help to end up 
with reliable conclusions in terms of the Hanle effect. 

CLASP will be observing the solar atmosphere for ap¬ 
proximately five minutes during which the structure of 
the dynamical TR will change to some extent. In our 
numerical experiments, we have only considered an in¬ 
stantaneous state of the solar model atmosphere corre¬ 
sponding to a particular snapshot of the time-dependent 
simulation. Furthermore, we have assumed an infinitely 
narrow spectrograph’s slit and we have neglected any 
possible drift of the slit during the simulated observa¬ 
tion. Because of the fine structuring and dynamical be¬ 
havior of the solar chromosphere, CLASP cannot aim 
at determining the magnetic field at high spatial reso¬ 
lution. However, the lack of spatio-temporal resolution 
might actually help us to end up with observables use¬ 
ful to determine the mean field strength of the solar TR. 
Therefore, the standard deviations shown in Figs. 14 and 
15 should be understood as being the most pessimistic 
estimates and the real observation will be of better ac¬ 
curacy. 

7. CONCLUSIONS 

The outer solar atmosphere is a highly structured and 
dynamic plasma, much more complex than a uniform, 
plane-parallel ID configuration. Very likely, the height 
of the TR varies from one place in the atmosphere to an¬ 
other, delineating a highly crumpled surface (e.g., such as 
in the top left panel of Fig. 3). Strong spectral lines origi¬ 
nate in this type of spatially complex plasmas, where the 
horizontal atmospheric inhomogeneities break the axial 
symmetry of the radiation field. This symmetry break¬ 
ing produces changes in the scattering polarization sig¬ 
nals, which may compete with those caused by the Hanle 
effect. Therefore, the interpretation of spectropolarimet- 
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ric observations in chromospheric and TR lines requires 
solving a complex radiative transfer problem, namely the 
generation of spectral line polarization in realistic 3D 
models of the solar atmosphere. Likewise, since the in¬ 
tensity profiles are practically insensitive to the magnetic 
field vector, the true way to validate or refute the models 
themselves is by confronting synthesized and measured 
Stokes profiles, because the magnetic field information is 
encoded in the spectral line polarization. 

In this paper we have extended to the full 3D case 
previous radiative transfer investigations of the scatter¬ 
ing polarization in the hydrogen Lya line, which were 
carried out using ID (Trujillo Bueno et al. 2011) and 
2D (Stepan et al. 2012) models of the solar atmosphere. 
The 3D model of the extended solar atmosphere cho¬ 
sen for this investigation is representative of the physical 
conditions of an enhanced network region, with magnetic 
field lines that reach chromospheric and coronal heights 
(see Fig. 1). Using the radiative transfer code PORTA 
developed by Stepan & Trujillo Bueno (2013) we solved 
the full 3D radiative transfer problem of the generation 
and transfer of scattering polarization in the hydrogen 
Lya line, and carried out numerical experiments to un¬ 
derstand its magnetic sensitivity due to the Hanle effect. 

The first noteworthy point is that contained in Fig. 3, 
which shows the radiation field tensors that quantify the 
symmetry properties of the Lya radiation. The frac¬ 
tional anisotropy of the Ly a radiation (see the top right 
panel) is fully negligible throughout the whole photo¬ 
sphere and chromosphere, but it suddenly becomes sig¬ 
nificant (of the order of 10 %) and predominantly nega¬ 
tive right at the TR.^^ The other panels of Fig. 3 show 
that the radiation field tensors that quantify the breaking 
of the axial symmetry of the Ly a radiation are also neg¬ 
ligible below the TR; i.e., in the regions where the optical 
depth at the line center is larger than unity. As the op¬ 
tical depth at the Ly a line center becomes smaller than 
unity at the TR, such radiation field tensors become sig¬ 
nificant. As a result, the emergent Lya radiation shows 
conspicuous scattering polarization signals, with sizable 
Q/I and U/I line-center amplitudes all over the solar 
disk, even in the unmagnetized reference case. 

At full spatial resolution the most likely fractional po¬ 
larization amplitude is 1.5% for simulated observations 
at /r = 0.4 (close to the limb scattering geometry), while 
it is 0.7% at /i = 1 (forward scattering geometry). The 
polarization signals drop below the 1% level when the 
spatial resolution is degraded. For example, for a square 
resolution element of 3 arc seconds the most likely for¬ 
ward scattering polarization signal decreases by a factor 
two. In the limit of no spatial resolution we find that the 
CLV of the line-center Q/I signal is qualitatively similar 
to that found in semi-empirical ID models of the so¬ 
lar atmosphere, with the largest polarization amplitudes 
around /i = 0.4 and with at /x = 1. 

Although the average magnetic field strength of the 
model’s TR is only {B) ^ 15 G, the Hanle effect in the 
hydrogen Lya line (whose critical Hanle field is 53 G) 


produces a significant depolarization (see Fig. 13). For 
stronger magnetic fields the ensuing depolarization is 
larger, until it reaches the Hanle saturation limit for 
{B) > 250 G. As shown in Fig. 15, for the radial slit case 
the {Q)/{I) line-center amplitudes are slightly smaller 
than in the tangential slit case, but they are equally 
measurable and show a similarly significant Hanle de¬ 
polarization. In Section 6.4 we give our arguments in 
favor of the radial slit option when observing the Lya 
scattering polarization with the GLASP sounding rocket 
experiment. 

Another important conclusion is that the Hanle effect 
in forward scattering geometry (i.e., when observing the 
solar disk center) does not necessarily create linear po¬ 
larization in the presence of an inclined magnetic field, 
as it happens in the idealized case of ID model atmo¬ 
spheres. The fact that in a 3D model of the solar atmo¬ 
sphere the radiation field that illuminates each point of 
the TR plasma does not have axial symmetry with re¬ 
spect to the local vertical direction implies that even in 
the absence of a magnetic field forward-scattering pro¬ 
cesses can produce significant linear polarization signals 
in hy a. Interestingly, we have found and understood 
that in such a spatially-complex plasma the Hanle effect 
of an inclined magnetic field tends instead to destroy 
the forward-scattering polarization, as it happens when 
observing close to the edge of the solar disk where one 
reaches the case of 90° scattering geometry. This and the 
previously summarized results are important for a reli¬ 
able interpretation of future observations of the Ly a po¬ 
larization produced by optically pumped hydrogen atoms 
in the chromosphere-corona transition region of the Sun. 


Financial support by the Grant Agency of the Gzech 
Republic through grant P209/12/P741 and project 
RVO:67985815, as well as by the Spanish Ministry 
of Economy and Gompetitiveness through projects 
AYA2010-18029 (Solar Magnetism and Astrophysi- 
cal Spectropolarimetry) and CONSOLIDER INGENIO 
GSD2009-00038 (Molecular Astrophysics: The Herschel 
and Alma Era) is gratefully acknowledged. This re¬ 
search was also supported by the Research Gouncil of 
Norway through the grant “Solar Atmospheric Mod¬ 
elling” and by the European Research Gouncil under 
the European Union’s Seventh Eramework Programme 
(EP7/2007-2013) / ERG Grant agreement Nr. 291058. 
The radiative transfer simulations presented in this paper 
were carried out with the MareNostrum supercomputer 
of the Barcelona Supercomputing Center (National Su¬ 
percomputing Center, Barcelona, Spain). We gratefully 
acknowledge the computing grants, technical expertise 
and assistance provided by the Barcelona Supercomput¬ 
ing Center. We are also grateful to the European Union 
COST action MP1104 (Polarization as a Tool to Study 
the Solar System and Beyond) for financing short-term 
scientific missions at the lAC that facilitated the devel¬ 
opment of this research project. 


REFERENCES 

Belluzzi, L., Trujillo Bueno, J., & Stepan, J. 2012, ApJ, 755, L2 Carlsson, M., Hansteen, V., Gudiksen, B. V., Leenaarts, J., & De 

Carlin Ramirez, E., Asensio Ramos, A., & Trujillo Bueno, J. 2013, Pontieu, B. 2014, in preparation 

ApJ, 764, 40 Curdt, W., Tian, H., Teriaca, L., Schiihle, U., & Lemaire, P. 2008, 

A&A, 492, L9 


Scattering polarization of Ly a in a 3D model of the chromosphere-corona transition region 


15 


De Pontieu, B., et al. 2014, Sol. Phys., 289, 2733 

Fontenla, J. M., Avrett, E. H., & Looser, R. 1993, ApJ, 406, 319 

Gudiksen, B. V., Carlsson, M., Hansteen, V. H., Hayek, W., 
Leenaarts, J., & Martmez-Sykora, J. 2011, A&A, 531, A154 

Ishikawa, R., Asensio Ramos, A., Belluzzi, L., Manso Sainz, R., 
Stepan, J., Trujillo Bueno, J., Goto, M., Tsuneta, S. 2014, ApJ, 
787, 159 

Kano, K., et al. 2012, Proc. SPIE, Vol. 8443, article id. 84434F, 17 
PP 

Kobayashi, K., et al. 2012, ASP Gonf. Ser., 456, 233 

Kneer, F. 1981, A&A, 93, 387 

Kubo, M., et al. 2014, in Solar Polarization 7, ASP Gonf. Ser., 489, 
307 

Landi DegPInnocenti, E. &: Landolfi, M. 2004, Polarization in 
Spectral Lines^ Kluwer Academic Publishers. 

Leenaarts, J., Garlsson, M., Hansteen, V., & Rutten, R. J. 2007, 
A&A, 473, 625 


Leenaarts, J., Garlsson, M., & Rouppe van der Voort, L. 2012, ApJ, 
749, 136 

Manso Sainz, R., & Trujillo Bueno, J. 2011, ApJ, 743, 12 
Roussel-Dupre, D. 1982, ApJ, 256, 284 
Stepan, J., &: Trujillo Bueno, J. 2011, ApJ, 732, 80 
Stepan, J., Trujillo Bueno, J., Garlsson, M., &: Leenaarts, J. 2012, 
ApJ, 758, L43 

Stepan, J., & Trujillo Bueno, J. 2013, A&A, 557, A143 
Truiillo Bueno, J. 2014, in Solar Polarization 7, ASP Gonf. Ser., 
489, 137 

Trujillo Bueno, J., Landi Degl’Innocenti, E. 1997, ApJ, 482, L183 
Trujillo Bueno, J., Stepan, J., & Gasini, R. 2011, ApJ, 738, Lll 
Trujillo Bueno, J., Stepan, J., & Belluzzi, L. 2012, ApJ, 746, L9 
Vourlidas, A., et al. 2010, Sol. Phys., 261, 53 




